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SELF-PERPETUATING SPIRAL ARMS IN DISK GALAXIES 

Elena D'Onghia 1 ' 2 ' 3 , Mark Vogelsberger 1 , Lars Hernquist 1 

ABSTRACT 

The causes of spiral structure in galaxies remain uncertain. Leaving aside the grand bisymmetric spirals with 
their own well-known complications, here we consider the possibility that multi-armed spiral features originate 
from density inhomogeneities orbiting within disks. Using high-resolution N-body simulations, we follow the 
motions of stars under the influence of gravity, and show that mass concentrations with properties similar to 
those of giant molecular clouds can induce the development of spiral arms through a process termed swing 
amplification. However, unlike in earlier work, we demonstrate that the eventual response of the disk can 
be highly non-linear, significantly modifying the formation and longevity of the resulting patterns. Contrary 
to expectations, ragged spiral structures can thus survive at least in a statistical sense long after the original 
perturbing influence has been removed. 

Subject headings: galaxies: kinematics and dynamics, galaxies: spiral, methods: numerical 



1. INTRODUCTION 

Seventy percent of galaxies in the nearby Universe are char- 
acterized by a disk with prominent spiral arms, but our under- 
standing of the origin of these patterns is incomplete, even af- 



ter decades of theoretical study ([Toom re 1977; Ath anassoula 
T984} |Binney & Trema ine 2008[ |Sellwood||2011| ). Several 



ideas have been proposed to explain the formation of spiral 
arms. One model posits that these features are large-scale 
density waves continuing to propagate in a differentially ro- 
tating disk. In particular, this theory argues that the matter in 
the galaxy (stars and gas) can maintain a density wave through 
gravitational interactions even in the presence of shear. This 
density wave remains at least quasi- stationary in a frame of 
reference rotating around the center of the galaxy at a fixed 
angular speed, identified wit h the pattern speed of the spirals, 
and c overs the entire disk ( [Lin & Shu||1964[ |Bertin & Lin| 
[T996] ). 

An alternative theory proposes that spiral arms are stochas- 
tically produced by loca l gravitational amplification in a dif- 
ferentially rotating disk ( [Goldreich & Lynden-Bell|1965[ |Ju-| 
|lian & Toomre 1966]). The mechanism behind this process is 
known as swing amplification and it can be seeded either by 
preexisting leading waves or else by the response of a disk 
to the presence of a corotating overdensity, such as a giant 
molecular cloud. This dynamical response takes the form of 
wakelets in the surrounding medium, each amplified by its 
own self-gravity through the swinging of leading features into 
trailing ones owing to the shear. 

According to this second theory, spiral arms would fade 
away in one or two galactic years if the dr iving perturbations 
were removed (Toomre & Kalnajs 1991), in contrast to the 
quasi-steady nature of the arms in the model proposed by Lin 
and Shu. Thus, a continuous source of perturbations would be 
required for these fluctuating spiral patterns to be maintained 
throughout the lifetime of a galaxy. Indeed, by mimicking 
the ef fects of dissipative infall of gas, Sellwood & Carlberg 
1984j > showed that the addition of fresh particles o n cir cular 



orbits could cause such spiral patterns to recur, and Carlberg 
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& Freedman ( 1985 ) demonstrated that almost any mechanism 
of dynamical cooling can maintain spiral activity of a similar 
kind. 

However, the pioneering work on density waves in the 
1960s was based on analytic theory and invoked various as- 
sumptions in order for solutions to be found. For example, the 
swing amplification analysis involved linear approximations 
to the equations of motion. The theoretical emphasis since 
that time has been on identifying driving mechanisms that can 
sustain the wave in spite of the damping which would cause 
it to decay in this picture. Indeed, whereas observations show 
that spiral arms might be density waves, N-body experiments 
have not yielded long-lived spiral structures, as predicted by 
the stationary density wave theory. Simulations of cool, shear- 
ing disks always exhibit recurrent transient spiral activity and 
this situation has not changed over th e past several decades 
as computational power ha s increased fSellwood|2000 Fujii 
|et al.|2011t[Sellwood|2011| ). 

Some work showed that spiral patterns fade away in numer- 
ical simulations of stellar disks if the effects of gas dissipation 
are not included; the reason is that the disk becomes less re- 
sponsive as random motions rise owing to partic le scattering 
by the spiral activity and giant molecu lar clouds ( Sellwood & 
|Carlberg|1984[|Binney & Lacey|1988| ). Moreover, the debate 
about the longevity of the arms practically ceased two decades 
ago because the available computational power did not permit 
definitive tests of some of the predictions of the theories and 
also because observations at that time were not sufficiently de- 
tailed to discriminate between the two main competing views. 

In the past decade, some studies have argued that the contin- 
uous infall of substructures in the dark matter halos of galax- 
ies could induce spiral patterns in disks by generat ing local- 
ized disturbances that grow by swing amplification ( Gauthier 
|et al.||2006} |Purcell et al.||201l| . According to such simula- 
tions, the main agent producing transient features would be 
satellite passages through the inner part of a disk. Because 
the tidal effects of the satellites are generally small (D'Onghia 



et al.|2010b| ), this process is distinct from interactions thought 
to be responsible for grand-design spirals like M5 1 . However, 
there are indications that dark matter substructures orbiting 
in the inner regions of galaxy halos would be destroyed by 
dynamical processes such as disk shocking, and he nce would 
not be able t o seed the formation of spiral structure ( [D'O nghia 
|et al.|2Q10a| ). 
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In what follows, we report on high resolution N-body sim- 
ulations of isolated disks in which we follow the motion of 
100 million star particles under the influence of gravity, to 
explore the idea that spiral arms are seeded by density inho- 
mogeneities orbiting within the disk. These inhomogeneities 
can be identified with fluctuations in the distribution of gas in 
the interstellar medium of galaxies, such as giant molecular 
clouds or, more speculatively, any massive star clusters em- 
bedded in the disk. Because early works already showed that 
the gravitational response of a stellar disk to an overdensity 
orbiting in its plane is similar a s the response of a gaseous 
disk (see in particular Fig. 7 of Toomre] ( |1981| )), we perform 
our analysis only on stellar disks, for simplicity. 

Our work focuses on understanding the origin of spiral 
arms and differs in many respects from previous studies. The 
methodology we employ here has the potential to discern the 
physical processes occurring in stellar disks at a higher level 
of detail than in previous simulations of isolated galaxies. 

As we describe below, our results indicate that the response 
of a disk to local perturbations is highly non-linear and time- 
variable on galactic scales and so is not fully captured by the 
linear approximations invoked in swing amplification theory 
or the classic work on quasi- steady density waves. 

In what follows, we describe our methodology in §2. In §3 
we present our results by showing examples in which a disk 
is perturbed by inhomogeneities orbiting in its plane, and we 
even vary the lifetimes of the perturbers. Results anticipated 
by linear theory are described in §4, whereas §5 is devoted to 
illustrating the non-linear effects that arise in our simulations. 
Finally, we discuss the astrophysical relevance of our work in 
§6. 

2. METHODOLOGY 

Our simulations were carried out with the paralle l TreePM 
code GADGET- 3 (last described in |Springel|20 05). We only 
employ the tree-based gravity solver coupled to a static exter- 
nal potential to solve for the evolution of collisionless parti- 
cles. 

P airwise particle interactions are softened with a spline ker- 
nel ([Hernquist & Katz 1989) on scales h s , so that forces 
are strictly Newtonian for particles separated by more than 
2h s . The resulting force is roughly equivalent to traditional 
Plummer- softening with scale length e « h s /2.S. For our 
simulations the Plummer-equivalent gravitational softening 
length set to e = 5pc for stellar particles. 

2.1. Setting up initial conditions 

The galaxies in our study consist of dark matter halos and 
rotationally supported disks of stars. The parameters describ- 
ing each component are independent and the models are con- 
structed in a manner similar to the approach described in pre- 
vious works ( |Hernquist||1993[ |Springel|2000[ jSpringel et aL] 
|2005] ). 

2.2. Dark Halo 
We model the dark matter mass distribution with a IHern-l 



quistlp990) profile: 
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length and M dm is the total halo mass here set to 9.5x1 11 
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In the past, models of disk galaxies run in isolation and used 
to study the properties of spiral arms employed only a few 
million particles to sample both the stellar disk and the dark 
matter halo. In such experiments, randomly-placed particles 
produce fluctuations in the halo potential. Even if the disk 
is initially featureless, the Poisson noise owing to such dis- 
cretization of the mass in N-body experiments is inevitably 
swing amplified, producing trai l ing multi-armed spiral pat- 
terns in the disk (Toomre 1977; Fujii et al. 2011; Sellwood 
[20121 ). 

In order to suppress the development of artificial features 
in all the N-body experiments that follow, we set up a live 
disk of stars embedded in a rigid dark matter potential. We 
employ simulations with a sufficiently large number of parti- 
cles in the disk, i.e. 100 million, so that the disks are essen- 
tially featureless when evolved without any perturbers acting 
on them. These simulations serve as "controls," making it 
possible to identify the response of the disk to imposed per- 
turbations. In this manner, we will be able to separate the 
sources responsible for exciting features in the disk from the 
stars which react to the perturbations, unlike previous experi- 
ments in which the stars themselves acted as perturbers, com- 
pli cating the interp retation of the experiments, as emphasized 
by |Toomre| ( ^990| ). 



2.3. Stellar Disk 

We model the stellar disk in the initial conditions with an 
exponential surface density profile of scale length R s : 



E*(i?) 



2irR 2 



exp(-R/R s ), 



(2) 



so that the stellar disk mass is M* = m^M tot , where is 
the disk fraction of the total mass M tot of the galaxy. For an 
extremely thin disk the circular velocity of the galaxy would 
be: 
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Here, G is the gravitational constant, y = R/(2R S ), and I n 
and K n are modified Bessel functions. We specify the vertical 
mass distribution of the stars in the disk by giving it the profile 
of an isothermal sheet with a radially constant scale height zo . 
The 3D stellar density in the disk is hence given by: 
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We treat z$ as a free parameter that is set by the vertical ve- 
locity dispersion of the stars in the disk and fix the veloc- 
ity distribution of the stars such that this scale height is self- 
consistently mainta ined in the full 3D potential of the galaxy 
| |Springel et al.|2005) . We adopt z = 0.1R S . 

Fig. [1J shows the rotation curve for the fiducial galaxy 
model used in the simulations (left panel). The total disk mass 
of the galaxy is 1.9xlO lo M . The total disk fraction is fixed 
to rrid = 0.02 and the scale length of the disk, R s , is 3.13 
kpc. The disk is set up so that it is stable, as measured by the 
Q parameter, defined for infinitely thin disks as: 

Q= 3.36GE' (5) 
where cr R is the radial velocity dispersion, k is the epicycle 
frequency, and E is the stellar surface density. The proper- 
ties of the disk are chosen so that the Q parameter is initially 
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Fig. 1 . — Left Parce/.Rotation curve for fiducial disk galaxy model (left panel). Right Panel. Toomre parameter Q as a function of radius for the fiducial model. 
The radius is normalized to the value of the disk scale length of 3.13 kpc in both panels. 



larger than 1 at all radii and has a minimum value of 1.3 (right 
panel of Fig[T]), implying that t he disk is stable to axisymmet- 
ric instabilities (Toomre 1964]). 

The strength of spiral features is estimated by dividing the 
face-on disk into concentric annuli which are further divided 
into azimuthal bins. The mass from the star particles is then 
assigned to these bins, making it possible to compute the stel- 
lar surface density E(i2, 0). Fourier components of the sur- 
face density are then calculated. We use the Fourier transform 
analysis to estimate the residuals of the surface density com- 
puted by subtracting the azimuthally averaged surface density 
from the surface density distribution and normalizing accord- 
ing to Res = (E - E avg )/E avg . 

3. RESULTS 

3.1. Disk Stability Test 

We first perform an N-body experiment in which a live disk 
of 100 million stars as described above is embedded in a rigid 
potential for the dark halo. We run the fiducial case with no 
perturbers for two-three galactic years and verify that the disk 
does not develop prominent spiral features. The outcome of 
this experiment is shown in Fig. [2] where a time evolution of 
the disk is displayed. 

3.2. Dynamical Response of the Disk to Perturbers 

Having established that relatively low mass disks run with 
sufficiently large numbers of particles do not develop promi- 
nent spiral patterns in the absence of any perturbing influence, 
we then evolve the same galaxy models, but add 1, 000 grav- 
itationally softened particles (each with a mass of the order 
of a typical giant molecular cloud: M = 9.5xl0 5 M ) dis- 
tributed with the same exponential law as the disk and as- 
sumed to be corotating on circular orbits with the disk stars. 
We note that we use a 20 times larger gravitational softening 
length for these molecular cloud particles due to their larger 
mass compared to the stellar particles in the disk. These per- 
turbers are kept moving on circular orbits in order to avoid 
any clumping that could facilitate the formation of overdense 
stellar regions. 

Each softened particle perturbs the motions of neighboring 
stars and drives a strong local response in the disk. The col- 



lective influence of all the perturbers leads to the development 
of spiral structures as shown in Fig. [3] These features extend 
well beyond two scale lengths of the disk. At a distance of 
five scale lengths these segments tend to be more isolated and 
weaker. 

In particular, Fig. [3] displays the time evolution of the disk 
for approximately two-three galactic years. Note that the dy- 
namical response is rapid and already after almost half of a 
galactic year, multi-armed patterns appear. Fig. [4] (top pan- 
els) shows the disk and the 1,000 molecular clouds distributed 
within it after approximately one galactic year. The disk is 
displayed face-on in Cartesian coordinates (left panels) and 
surface density residuals are shown in polar coordinates (right 
panels). In this case, interior to one disk scale length the disk 
is characterized by three or four spiral arms. At two scale 
lengths from the galaxy center the number of spiral arms is 
approximately seven with amplitudes reaching values of 10- 
15% higher than the stellar background. 

We repeat the same experiment with the giant molecular 
clouds counter-rotating relative to the disk stars. As shown in 
Fig. [4] (bottom panels), the disk response after approximately 
one galactic year is now much weaker than for corotating per- 
turbers displayed at the same time (top panels). 

3.3. Finite Lifetime Giant Molecular Clouds 

Our previous calculations assumed an infinite lifetime for 
the giant molecular clouds. Although the lifetimes of these 
systems are much debated, giant molecular clouds in galaxies 
are likely short-lived. Furthermore, giant molecular clouds 
in real galaxies are preferentially found in spiral arms, and 
may be created there by the converging flow of the interstellar 
medium due to its motion through the spiral potential. The 
clouds disperse after they pass through the arm, perhaps in 
part because of the diverging flow as the gas emerges from 
the arm, although other factors such as of the energy input 
from young stars are important also. This is certainly a too 
complicated picture to be simplified, however here we only 
want to measure the dynamical response of the disk to some 
perturbers with different lifetimes. 

In the following, we show the outcome of runs performed 
with a finite lifetime for the giant molecular clouds. In partic- 
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Fig. 2. — Time sequence for the evolution of a live disk of stars embedded in a dark Milky Way-sized halo run unperturbed in isolation and displayed after a 
few galactic years. The disk is displayed face-on in Cartesian coordinates. 



ular, we show two examples: a first run adopting the fiducial 
disk of stars embedded in a rigid halo potential with 1,000 
giant molecular clouds having a lifetime of 2.5xl0 6 yrs and 
a second run with cloud lifetimes of 2.5xlO T yrs. The fi- 
nite lifetime is directly implemented in our simulation code. 
Specifically, we assign a time counter to each giant molecular 
cloud particle. Once this counter exceeds the predefined life- 
time, this particle is converted to a regular stellar particle by 
changing its mass and gravitational softening length. We then 
randomly select another stellar particle in the disk and turn 
it into a giant molecular cloud particle also by changing the 
mass and gravitational softening length accordingly. To avoid 
the sudden destruction of all giant molecular cloud particles 
at the same time, we generate the initial set of giant molecular 
cloud particles with a uniform age distribution in the range of 
the selected lifetime. 

The top panels of Fig. [5] show the outcome of the run with 
the giant molecular clouds having a lifetime of 2.5xl0 6 yrs. 
This lifetime is too short to trigger any resonance between the 
shear flow and epicyclic vibration, hence the forced response 
that leads to the formation of the arms does not occur. The op- 
posite situation is shown in Fig. [5] (bottom panels) where the 
giant molecular clouds are assumed to live for 2.5xl0 7 yrs. 
Note that this lifetime is sufficient to trigger the formation of 
multi-armed spiral patterns. The orbital period for a typical 
giant molecular cloud in a Milky- Way sized galaxy is w 10 8 
yrs, thus our calculations show that a lifetime of a quarter of 
the orbital time is sufficient to trigger a strong response in a 
disk and the subsequent development of spiral arms. The ar- 
gument can be now inverted: in order to form realistic spiral 
arms in disk galaxies, the giant molecular clouds have to sur- 



vive for at least 10 7 yrs, which is in accord with observational 
estimates. 

4. RESULTS ANTICIPATED BY LINEAR THEORY 

4. 1 . Material Arms or Density Waves ? 

A long-standing controversy over the nature of spiral arms is 
whether they correspond to density enhancements in the back- 
ground stellar distribution (density waves), or are made up of 
stars that always remain in the arm and are just more con- 
centrated than the stars outside the arm (material arms). In 
the early studies, the arms were assumed to be density waves, 
because if they were material they would quickly wind up as 
the galaxy rotates. Thus, both the swing amplification and the 
static density wave theories argue that the arms are overdense 
regions of the disk moving around at a different speed rela- 
tive to the stars themselves. Stars thus continuously move in 
and out of the spiral arms. However, recent investigations us- 
ing numerical simulations of stellar disks have challenged this 
claim and and argu e that the arms might be material structures 
( [Grand et al.|2012| ). 

To investigate this in the context of our simulations, we 
identify a patch of stars along the arm in the stellar disk after 
the arms are fully developed (after 100 Myrs) as displayed in 
top panels of Fig. |6| where the patch is colored in black. Then, 
we follow the positions of the stars originally in the patch for- 
ward in time and display the outcome after two galactic years 
(bottom panels). We note that the stars initially in the patch 
spread out, confirming that the spiral patterns in our simu- 
lations are density waves and not material structures. This 
is shown in polar coordinates in the bottom panel of Fig. [6] 
where it is clear that the patch is being sheared out by differ- 




Fig. 3. — Time sequence of an N-body experiment in which a self-gravitating disk of 100 million stars seen face-on embedded in a dark Milky Way-sized 
halo is run with 1,000 giant molecular clouds (represented with the black dots) which are randomly distributed within the disk and are assumed to be corotating 
on circular orbits with the disk stars. The giant molecular clouds act as perturbers and the live disk dynamically responds to the presence of these perturbers by 
developing features which resemble multi-armed structures in galaxies. Each panel displays a region 30 kpc on a side at the times indicated. 

ential rotation and the pitch angle of the patch differs signifi- 
cantly from that of the spiral features. 



4.2. Forced-Wake Theory 

The different response of the disk to the direction of the 
motion of the perturbers within the disk, as shown in Fig. [4] 
suggests that the physical process behind spiral arm formation 
in our simulated galaxies is initially, at least, wake-making, 
which operates through a combination of three ingredients: 
the shearing fl ow, small-scale epicyc lic motions, and the disk 
self-gravity ( Julia n~& Toomre| [T966). We briefly review the 
main concept. 

Consider stars drifting in the flow at slightly larger radii 
than a perturber which is corotating with the disk. The drift- 
ing stars will feel a small force owing to the presence of the 
perturber, exciting epicyclic motion in the stellar orbits. At 
the same time, the inwards and outwards motions of the stars 
lead to density enhancements and deficits relative to the stellar 
background. These enhancements can swing from leading to 
trailing owing to shear in the disk. The sense of this shearing 
motion is identical to that of the retrograde epicyclic motion 
of individual stars. Because of this match it is possible for a 
swinging pattern to resonate with the epicyclic stellar motions 
as the shear causes it to become a trailing pattern. This tempo- 
rary match allows the feature to be amplifi ed by self-gravity 
as the stars linger within the perturbation ( Julian & Toomre 

[1955) . 

4.3. The Critical Role of Self- gravity as Amplifier 



Juli an & Toomre| ( |1966| ) showed that swing amplification 
results from a cooperative effect between shear, epicyclic mo- 
tions, and self-gravity. In order to verify the importance of 
self-gravity in the outcomes of our simulations, we performed 
an N-body simulation of a disk of stars treated as test particles 
embedded in a dark halo along with 1 ,000 perturbers moving 
on circular orbits within the disk (for details of the set-up, see 
Appendix). 

Fig. [7] shows a time sequence from a simulation in which 
the disk is perturbed by density inhomogeneities, but the self- 
gravity of the disk is ignored. The influence of the disk self- 
gravity is clear. In the absence of self-gravity, wakes still form 
around the perturbers, but the wakes have different shapes 
and are weaker than when self-gravity is included, as demon- 
strated in Fig. [7] A harmonic analysis shows that the magni- 
tude of each Fourier coefficient is greatly reduced if disk self- 
gravity is ign ored, as expected on the basis of similar tests 
performed by Juli an & Toomre | (|1 9 66|). (In p articular, com- 
pare Figs. 7 and 12 in | Julian & Toomre ( |1966| ).) Thus, self- 
gravity acts as a communicator between the stars, amplifying 
wakes around each perturber and allowing different wakes to 
connect to produce long-range patterns, as in Fig. [3] 

4.4. Collective Amplification Regime 

Swing amplification is the linear gravitational response of 
a disk to an overdensity orbiting in its plane. In a stable 
disk, pressure-like or mixing effects of random motions avoid 
clumping and normally disrupt an overdense region before 
it has time to collapse. Because of the similarity between 
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Fig. 4. — Top Panels. A live disk of stars embedded in a dark Milky Way-sized halo run in isolation with 1,000 giant molecular clouds distributed within the 
disk and assumed to be corotating on circular orbits with the disk. The disk is displayed approximately after one galactic year face-on in Cartesian coordinates 
(left) and surface density residuals are shown in polar coordinates (right). Bottom Panels. The same experiment is displayed at the same time with the 1,000 
perturbers assumed to be counter-rotating on circular orbits with respect to the disk. 

shear and epicyclic vibrations however, individual stars re- 
main within shearing density enhancements for a substantial 
fraction of the epicyclic frequency k, which is given by 
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an 2 

dR 
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where ft = V c /R is the angular frequency, R is the radius, 
and V c is the total circular veloci ty. For the ide alized case of 
a Mestel disk, with V c constant, Toomre (1981) found that the 
gain of the swing amplifier depends on the value of Q and the 
ratio of the wavelength to a critical value X — A/ \ cr it where 
A C rit = 47rGEft -2 is the shortest wavelength stabilized by ro- 
tation alone. For this idealized situation, the amplification 
factor peaks at X ~ 1.5, reaching values in excess of 100 
for cold disks. Swing amplification is negligible for X > 3. 



Indeed perturbations with such long wavelengths are stabi- 
lized by the angular momentum of the stars. Note that in 
N-body experiments, particle noise creates a spectrum of per- 
turbations of all sizes and shapes, including both leading and 
trailing spirals. As the differential rotation of the disk shears 
leading spirals into trailing configurations, swing amplifica- 
tion boosts the amplitudes of those which have wavelengths 
of order of A ~ 1.5A cr ^ thereby creating a multi-armed pat- 
tern of trailing spirals with a characteristic spacing that ex- 
plains the results illustrated in Fig. [3] (see Appendix for an 
expression giving the amplification parameter for a general 

disk galaxy). 

A naive expectation from | Julian~& Toomre (1966) is that 
each inhomogeneity in the disk should generate a wake, form- 
ing an initial leading feature which is amplified owing to the 
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Fig. 5. — Left Panels. The outcome of the N-body experiment where a live disk of stars embedded in a dark Milky Way-sized halo run with 1,000 giant 
molecular clouds with a finite lifetime of 2.5xl0 6 yrs (top) and 2.5xl0 7 yrs (bottom). Right Panels. The surface density residuals are shown in polar coordinates. 



differential rotation of the disk to become a stronger trailing 
wave. However, in our simulations the number of arms is 
not determined by the number of perturbers. 4 Instead, each 
perturber produces a segment of an arm and then these seg- 
ments are joined at kinks. Thus, there is a collective process 
that occurs between the individual wakes that connects them 
together to produce a pattern that resembles the features in 
multi-armed spiral galaxies. 

According to swing amplification theory, the most strongly 
amplified features are those which have wavelength ~ X C rit- 
If in some region the perturbers are packed together more 
tightly than this, they will effectively act as an individual per- 
turber producing a single response with an extent ~ X C rit- In 
regions where the perturbers are separated by distances much 

4 It should be noticed that in linear theory, there is fundamentally no differ- 
ence, whether or not the perturbers are distinct from the supporting medium 
{Toomre|1990} . 



larger than A cr ^, they are effectively isolated, and we would 
expect to see individual, disconnected wakes around each per- 
turber. However, in the intermediate regime, when the density 
of the perturbers is such that they can each produce wakes 
that interact with one another, the wakes will connect owing 
to self-gravity, joining together to produce patterns that can 
cover an entire disk. This has the effect of selecting out par- 
ticular global structures that extend throughout a disk and the 
resulting number of arms is governed by the structural prop- 
erties of a galaxy through A cr ^. 

4.5. Determination of the Number of Arms 

So far we have presented models of disks of such mass 
that when they are perturbed they develop features charac- 
teristic of multi-armed spiral galaxies. We have verified that 
the number of arms does not depend on the number of per- 
turbers. However, the response does depend in detail on the 
mass distribution of the galaxy and the extent to which the 
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Fig. 6. — Top Panels. A patch of stars is identified in the overdense region of the arm and colored in black at the time when the arms are fully developed (100 
Myrs). Bottom Panels. The stars initially identified in the patch are followed forward in time and displayed after less than two galactic years. 



disk is self-gravitating. 

In a set of experiments shown in Fig. [8] we altered the 
contribution of the disk to the total mass of the galaxy. We 
kept the total mass of the galaxy and the mass distribution 
of the dark halo as previously but we ran simulations where 
we increased the disk fraction, rrid, from 0.01 to 0.02, 0.03 
and 0.04. Once again, we first ran these models without per- 
turbers orbiting within the disk for one galactic year and ver- 
ified that the disks do not develop prominent patterns. We 
subsequently introduced in the initial conditions 1,000 giant 
molecular clouds corotating with the disk stars. By changing 
the disk mass these models were designed to vary the criti- 
cal length scale parameter X crit = Att 2 GTJj£ and lead to a 



different sp i ral morphology (see Fig .l in Carlberg & Freed- 
|man|fll985fr ; |Efstathiou et al.| ( 11982) ), from the multi-armed 
features typical of galaxies like NGC 7217 obtained for low- 



mass disks, to those with a few prominent arms as in M 101 or 
the Milky Way. As anticipated based on the above arguments, 
the results shown in Fig. [8] demonstrate that the spiral mor- 
phology in our simulated disks is indeed determined by the 
structural properties of the galaxy through A cr ^. 

5. RESULTS NOT ANTICIPATED BY LINEAR THEORY 

5.1. Self -perpetuating spiral arms 

Theories based on a linear approximation describing the 
local dynamic response, e.g. a disk patch, predict that the 
emerging spiral patterns should be a superposition of swing- 
ing wavelets, shearing with the general flow, but exhibiting 
transient growth as they swing from leading to trailing. If the 
perturbation responsible for exciting a wakelet is removed, 
the linear theory predic ts that the wakelet will d ecay, reveal- 
ing its transient nature ( jToomre & Kalnaj s"[T991| ). In order to 
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Fig. 7. — Time evolution of a disk perturbed by 1,000 inhomogeneities co-rotating with the stars, but in which the disk self-gravity is not included. 



test the validity of the linear theory and to determine if spi- 
ral structures arise collectively through a local response of the 
disk, we stacked the giant molecular clouds at a fixed radius in 
the disk, superposing their associated wakes. The outcome is 
surprising, and we find that the response is highly non-linear 
in the sense that the wakes that initially formed around the 
perturbers, e.g. the giant molecular clouds, depart from them 
after only about a tenth of a galactic year and then become 
new perturbers that continuously excite the disk. 

In fact, this is apparent in the top right panel of Fig. [4] 
where it can be seen that at a fixed radius between say 2 
and 4 disk scale-length, the wakelets have different shapes 
and cannot be superposed. Furthermore, the perturbers (the 
black filled circles) are associated with their each individual 
wakelets, as predicted by linear swing amplification theory, 
only at distances beyond 4 disk scale-lengths where the self- 
gravity of the disk is lower and the individual perturbers are 
well-separated. As described below, we find that between two 
and four disk scale-lengths the arms are reinforced owing to 
non-linear effects and the wakelets become sufficiently over- 
dense that they decouple from the original perturbers and then 
act as perturbers themselves. 

To demonstrate this unexpected outcome and its radical 
consequences, we remove the original perturbers from the 
disk by replacing the giant molecular clouds with an equal 
mass in stars, and continue to evolve the disk, as illustrated 
in Fig. [9] The patterns do not disappear, contrary to the ex- 
pectations of linear theory, and are instead long-lived. Be- 
cause of the self-gravity of the disk, the wakes imprinted by 
the original perturbers become mildly self-gravitating and are 
of sufficiently high density to serve as perturbers themselves, 



thereby making the spiral features self-perpetuating. We have 
verified that the ensuing spiral patterns survive for at least an- 
other eight galactic years. 

We emphasize that in our simulations, the spirals are self- 
perpetuating in a manner completely distinct from the classi- 
cal density wave theory where the patterns are static or slowly 
rising in amplitude with tim e and are rigidly rotating ( |Lin &| 
|Shu| 1964[ |B"ertin & Lin|1996] ). Indeed, in our simulations the 
patterns are not global but they appear as local segments con- 
necting together and fluctuating in amplitude with time. The 
evolution of the resulting spiral structure is characterized by a 
balance between shear and self-gravity. Locally, shear tends 
to stretch the waves, breaking the arms, whereas in regions 
where the self-gravity dominates, the pattern is overdense and 
generates the segments making up the arms. This local bal- 
ance between shear and self-gravity gives the visual impres- 
sion that the spiral structures are global patterns, but in reality 
they are only segments, produced by local under-dense and 
over-dense regions. 

An alternative explanation for our results involves the pos- 
sibility that such dynamical systems support actual coher- 
ent unstable modes, and that the features observed after the 
perturbers are removed are manifestations of these unstable 
modes. Thus the rapidly varying spirals seen in our simula- 
tions would result from the superposition of a few long-lived 
waves (for a detailed study, see Sellwood ( 2012])). 

We generated spectrograms to test this interpretation, which 
supported at some level this possibility. However, our exper- 
iments were not run for a sufficiently long time to be deci- 
sive on this point. We will explore this alternative explanation 
with tests to discriminate between the two interpretations in a 



10 



D'Onghia et al. 




Fig. 8. — Sequence of models where we varied the contribution of the disk to the total mass distribution in the rotation curve of the galaxy. We fix the total 
mass of the galaxy and the mass distribution of the dark halo, but increase the disk mass fraction from a value of 0.01 to the following values: 0.02, 0.03 and 
0.04. Increasing the disk fraction decreases the number of arms. Models are displayed after one galactic year. 



forthcoming paper. 

Our results naturally predict that the pattern speed of spi- 
ral arms in late-type spiral galaxies should not be constant 
but should vary with distance from the center of the galaxy. 
This has already been claimed for patterns speeds estimated 
for N GC 5164 and NGC 628 JFoyle et al.|2010[ 



2009) and M101 ( [Foyle et al.|2011| ). HoweverTw] 



Meidt et al.l 



hereas these 

works conclude that this supports the interpretation that spiral 
arms are short-lived, our simulations show that the features 
can be long-lived as various local segments break and then 
later reconnect with other pieces of spiral arms, resulting in 
long-lived structures. 

Note that we have focused on the idealized case of a rather 
low disk mass. These systems are stable over long periods of 



time and do not develop bars or arms if unperturbed. If we 
consider disks that are more strongly self-gravitating, by, for 
example, increasing the ratio of the mass of the disk to the 
total mass of the galaxy, = 0.04, the simulations favor the 
formation of four spiral arms as displayed at the top left panel 
of Fig. 8] In that event, X crit is larger, so that longer wave- 
length disturbances are most strongly amplified. By examin- 
ing the rotation curves of our galaxy models, we also find that 
those with more massive disks have lower shear rates where 
the bulk of their mass is located, so individual wakes can re- 
main connected over longer distances without breaking apart 
owing to shear, yielding patterns with fewer global arms. In- 
deed the number of arms does not depend on the number of 
perturbers but on the structural properties of the galaxies in a 
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Fig. 9. — Time sequence for the evolution of spiral arms once the giant molecular clouds are removed from the disk. The left upper panel shows the stellar 
disk the time at which the perturbers, e.g. giant molecular clouds, are removed. We define this initial time as set to t=0. The subsequent time evolution of the 
spiral arms is shown in the other panels according to the initial time definition. Each panel displays a region 30 kpc on a side at the times indicated. 



way that by changing the rotation curve of the galaxy also the 
number of arms varies accordingly. 

In a further set of experiments, we ran simulations of disks 
after we systematically reduced the time at which the 1,000 
perturbers, e.g. giant molecular clouds, were removed and 
estimated the minimum time the perturbers need to corotate 
with the disk to trigger the formation of spiral arms that per- 
sists even after being removed. We find that if we remove the 
1,000 perturbers after 15 Myrs, this time is sufficient to trig- 
ger the process that leads later to the formation of long-lived, 
self-perpetuating patterns, similar to the estimates we made 
earlier for the required lifetimes of the original perturbers. 

5.2. Non-linear effects 

To demonstrate explicitly the possible importance of non- 
linear dynamics in spiral structure formation, the top panel 
of Fig. [10] shows the case of a galactic disk influenced by a 
single perturber with a mass of 10 7 M kept cycling at 2 disk 
scale length, corotating with the disk, and then being removed 
after one galactic year. 5 The image displays the disk when 
the perturber is still corotating with it. Initially, a wakelet is 
generated around the perturber in a manner consistent with 
the linear theory. However, non-linear effects are evident in 
the gravitational response of the disk to the perturbation and 
result in the development of additional underdense and over- 
dense regions that are clearly visible especially in the polar 
coordinate plot (right top panel of Fig. [TO]). These features, 



which act as perturbers, are responsible for keeping alive the 
spiral activity and are not expected based on the linear anal- 
ysis. These findings indicate that our simulations are probing 
an entirely new regime in which the linear theory is no longer 
fully applicable. 

Once the non-linear features begin to develop, we remove 
the perturber as shown at the bottom panel of Fig. 10 and 
find that the wakelets do not disappear, as predicted by the 
linear theory. Indeed, the overdense and underdense regions 
act themselves as new perturbers, exciting a further response 
in the disk, leading to a self-perpetuating pattern in which ad- 
ditional unexpected features such as holes and wrapping arms 
in the outer regions are produced within four galactic years, 
as shown in Fig. 
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5 The perturber is forced to cycle on a circular orbit so that it is not affected 
by dynamical friction. 



The expectations of linear theory are 
that the disk should relax after the perturbation is removed 
and that the spirals should gradually disappear, but we find 
instead that the patterns survive for many rotation periods. 

We note here that the response in the disks we simulate orig- 
inates at the co-rotation resonance between the perturbers and 
the surrounding stars. Thus, while we do not argue against 
the possi bility that dynamical ef fects such as "groove insta- 
bilities" ( Sellwood & Lin| [T989 ) can play an important role 
in producing spiral structure in actual galaxies, in our simu- 
lations, at least, the dominant dynamical effects seem to be 
associated with the response at the co-rotation resonance in 
the disk. 

Next, we performed the same experiment with a perturber 
mass reduced by a factor of ten, e.g. typical of a giant molec- 
ular cloud or a globular star cluster of mass 10 6 M kept cy- 
cling at two scale lengths of the disk and checked how long 
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Fig. 10. — Top Panels: A stellar disk responding to a single perturber with mass of 10 7 Mq, which is removed after one galactic year. Bottom Panels: Surface 
density is shown 4 galactic years after the perturber is removed confirming that the patterns survive for several rotation periods. The disk is displayed face-on in 
Cartesian coordinates (left panels) and surface density residuals are shown in polar coordinates (right panels). 



it takes for non-linear effects to develop. The outcome is il- 
lustrated in Fig. [TT] which shows that the linear approxima- 
tion breaks down after roughly 3 galactic years and holes and 
winding arms appear following an additional 2 galactic years 
once the perturber is removed. This suggests that while reduc- 
ing the mass of the perturber lengthens the time for non-linear 
effects to become important, they nevertheless still develop. 
Note that the amplitude of the arms (which are the overdense 
regions) and the depth of the holes (the underdense regions) 
are also reduced because they are sensitively dependent on the 
mass of the perturber. 

6. PREDICTIONS 

It seems clear from our results that a modest overdensity 
corotating with a galactic disk for a small fraction of a galactic 
year can trigger the formation of spiral structure that will take 



over and become self-perpetuating for at least several billion 
years. Intriguingly, observations of nearby disk galaxies sug- 
gest that when counter-rotating stellar or gas disks are present 
in galaxies there is no evidenc e of spiral structure ( [Bertola 
et al.|1999}|Coccato et al.|20TT] ), which is as expected, based 
on our simulations. 

Furthermore, it has been argued that spiral structures are 
transient features that need to be re-generated by continuously 
adding new perturbers in the form of gian t molecular clouds 
or clumpy gas that maintains ragged arms ( [Sellwood & Carl- 
berg|1984| ). Unlike simulation results of decades ago, which 



were compromised by limited numerical resolution, our ex- 
periments show that disks do not necessarily thicken because 
of the presence of spiral arms and giant molecular clouds, 
and that the spiral structure can persist for timescales simi- 
lar to the ages of actual galactic disks. Moreover, although 
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Fig. 11. — Top Panels: A stellar disk responding to a single perturber with mass of 10 6 Mq, which is forced to remain on a circular orbit (represented by 
the black circle). The image displays the disk after three orbital periods computed at a distance of two scale lengths of the disk where the first non-linear effects 
are visible. Bottom Panels: The patterns survive for several rotation periods and surface density is shown after 2 additional galactic years once the perturber is 
removed. The disk is displayed face-on in Cartesian coordinates (left panels) and surface density residuals are shown in polar coordinates (right panels). 

new gas accreted can excite young arms, there should be old 
stellar arms present in galaxies and fresh gas is not required 
to maintain spiral patterns over long timescales. This means 
that, in principle, early-type galaxies which are now gas-poor 
may still exhibit weak spiral features . IC 3328, an e arly-type 
dwarf galaxy with spir al structur e ( [Lisker &F uchs 2009]) or 
NGC1533 ( [DeGraaff et al.|2007) may be such cases. 

Ultimately, other predictions of our models include defi- 
nite quantitative statements regarding the pattern speeds and 
time-variability of actual spiral features in galaxies (Nelson, 
D'Onghia & Hernquist 2012, in prep.), and the kinematic and 
spatial response of the stars to these fluctuating features. Spi- 
ral waves in galaxy disks churn the stars and gas in a manner 
that largely preserves the overall angular momentum distribu- 
tion and leads to little increase in random motion. However, 



recent works ( Sellwood & Binney 2002 ) claim that changes in 
the angular momenta of individual stars are typically concen- 
trated around the co-rotation radius for an individual spiral 
wave, but given the supposed transient nature of the waves 
with a wide range of pattern speeds developing in rapid suc- 
cession, the entire disk should be affected. Our models indi- 
cate that spiral structures are not transient features, but even- 
tually locally break and subsequently reconnect as a balance 
between the shear that stirs the arm and the self-gravity that 
tends to reconnect the pattern. This calls for a better under- 
standing for how stars like the Sun can migrate if the nature 
and the evolution of t he waves is different ( Ro skar et al. 2008 ; 
|Minchevet al.|2011) 

We anticipate that there should be a correlation between the 
amplitude of the arms and the Q parameter that defines how 
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dynamically hot the disk is, in the sense that "colder" disks 
should, on average, display more prominent spirals. This 
follows from the fact that swing amplification, which is still 
central to our interpreta tion, is most vigorous in dynamically 
cool disks (flbo mre 1981). Ongoing and planned observations 
(Elmegre en et al.||2011| ) will enable tests of our modeling in 
the near future, which should shed further light on the still- 
mysterious nature of spiral structure. 

7. CONCLUSIONS 

In this paper, we have examined the possibility that spi- 
ral structure in disk galaxies originates through the dynamical 
response of self-gravitating, shearing disks to local density 
perturbations. While our analysis has much in common with 
earlier studies of spiral structure formation, it differs in the 
following respects. 

• Our N-body simulations show that spiral patterns can 
be described by density waves but the response of the 
disk to local perturbations can be highly non-linear and 
time- variable. Our work indicates limitations of two 
traditional theories that tried to explain arm formation: 
the t heory of arms as static density waves (|Lin & Shu 
1964) and linear forced- wake theory ( [Julian & Tbomre 



1966) which were based on simplifying assumptions. 
In particular, our results are inconsistent with the idea 
that spirals are static density waves because the arms in 
our simulations fluctuate in time. Moreover, our results 
are different from what is expected based on forced- 
wake theory because we show that the arms can be 
long-lived and are not necessarily transients even when 
the original perturbations responsible for driving the 
patterns are removed. 



• The spiral patterns in our simulations are not global as 
predicted by the classical static density wave theory but 
locally they appear to fluctuate in amplitude with time. 
It is the local balance between the shear that stirs waves, 



breaking the arms, and the self-gravity that makes up 
the arms that gives the visual impression that the pat- 
terns are global whereas they are actually segments pro- 
duced by local under-dense and over-dense regions, at 
least for relatively low mass disks. These under-dense 
and over-dense regions act as perturbers, keep alive 
the spiral activity and probe an entirely new regime in 
which the linear theory is no longer fully applicable. 

• The dynamical response in the disk originates at the co- 
rotation resonance between the perturbers and the sur- 
rounding stars. This implies that our features that wrap 
up and get replaced statistically have angular speeds 
that are not constant with location, but decrease with 
radius. The feature speed slightly differs from the sur- 
rounding stars, explaining why the arms do not wind up 
while being at the same time long-lived. 

• A modest overdensity corotating for at least 25 million 
years will excite a response in the disk that will lead to 
a self-perpetuating pattern owing to non-linear effects. 

We believe that there is already some observational support 
for our picture. In particular, our model predicts that there 
should exist red disk galaxies with no gas but that ma y still 
have arms. Gal axies such as this have been discovered (Mas- 
ters et al.|2010| ) and current theories for spiral arm formation 
do not explain their morphologies. In our model, the pattern 
speeds of the spiral features should vary with location within 
a disk, as has been claimed for som e nearby spirals ( |Foyle| 



|et al.|2010|[20TT||Meidt et al.|2009| ). We also predict that the 
number of arms in galaxies should depend on the structural 
properties of the galaxy and that there will be a correlation 
between the amplitude of the arms and how dynamically cold 
the disk is, in the sense that colder disks should display more 
prominent arms. These predictions are readily testable and 
relevant observational studies should provide strong tests of 
our model in the forthcoming years. 



APPENDIX 

STATIC POTENTIAL FOR A THIN EXPONENTIAL DISK 

To verify the importance of self-gravity in the formation of arms, we performed an N-body simulation of a disk of stars treated 
as test particles. In order to treat the stars as test particles, we compute the forces acting on them by introducing a static potential 
for a thin exponential disk of the form: 

a x = -cos6» — 
oR 

a y =-sme M 

GM * exp (-R/R h ) tanh (j\ (Al) 



where 9 is the angle phase and 



dr 



GM*z 



exp(-R/R h )(^-) lncosh(^) 



GM, 



[/ (y)tfi(y)-/i(y)tfo(y)] 



2R h 
GM*y . 



I {y)K[{y)- 



- I[(y)K (y) - h{y)K(y)\ (A2) 
The interactions between the perturbers and stars are computed directly, meaning that the stellar orbits are influenced by these 
gravitational interactions, as in the simulations with disk self-gravity presented above. 
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AMPLIFICATION PARAMETER FOR A GENERAL DISK GALAXY 
To interpret the results shown in §4.4 we consider the Toomre X -parameter for a galaxy like that in our simulations, using 

K 2 



X = 



^ 2 R 



(Bl) 



27rG£ra 

In general, the epicyclic frequency depends on the total angular frequency which is the sum of the angular frequencies of all the 
components of a galaxy, in particular the bulge, disk and dark halo, according to 



n 2 = n 2 D + n% 



n 2 H . 



For an exponential disk: 



2 _ GM D 



h{y)K (y) - I 1 (y)K 1 (y) 



(B2) 



(B3) 



where y = R/2Rf l , is the disk scale length and Md the disk total mass. For bulge and halos described by Hernquist models, 
the angular frequency of the bulge is 

2 _ GM B 

- (i?r^r (B4) 

where a 5 and Mb are the bulge scale length and mass, respectively and for the halo 

2 _ GM H 

n » - (RT^y (B5) 

where and Mb are the halo scale length and the mass, respectively. The amplification parameter for a general disk galaxy 
becomes: 



X 



m 



M B 2y + 3a b /R h 



M D (2y + a h /R h f 
M H 2y + 3a h /R h 



M D (2y + a h /R h y. 
- (shKo - 3/oifi + hK 2 - I 2 K^ 



Ay(I K - IiJTi)) 



(B6) 



where m characterizes the Fourier coefficient A m and indicates the number of arms. Note that substituting the values for the 
structural parameters of the galaxy model adopted here into eq. (B6) gives an estimate of the number of arms as ~ 7, in accord 
with our simulations. 
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